Geometric Phase Gates with Adiabatic Control in Electron Spin Resonance 
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High-fidelity quantum operations are a key requirement for fault-tolerant quantum information processing. In 
electron spin resonance, manipulation of the quantum spin is usually achieved with time-dependent microwave 
fields. In contrast to the conventional dynamic approach, adiabatic geometric phase operations are expected 
to be less sensitive to certain kinds of noise and field inhomogeneities. Here, we investigate such phase gates 
applied to electron spins both through simulations and experiments, showing that the adiabatic geometric phase 
gate is indeed inherently robust against inhomogeneity in the applied microwave field strength. While only little 
advantage is offered over error-correcting composite pulses for modest inhomogeneities < 10%, the adiabatic 
approach reveals its potential for situations where field inhomogeneities are unavoidably large. 



Precise coherent control of quantum systems is an essen- 
tial ingredient for many quantum technologies. In partic- 
ular, high-fidelity gate operations on quantum bits (qubits) 
are central to practical realisations of quantum information 
processing {]]]. The electron spin provides a quantum two- 
level system that is well-suited for the physical implementa- 
tion of a qubit. Dynamic control of electron spin states is 
commonly realized by applying microwave pulses in elec- 
tron spin resonance (ESR) 131 ■ Although single microwave 
pulses in conventional ESR spectrometers usually have non- 
negligible errors in amplitude and phase JH, high-fidelity 
single-qubit operations can often still be realized usingcare- 
fully designed composite pulses such as BB1 pulses 14j,|5[] and 
Knill pulses flU. 

A different approach to qubit operations involves geometric 
manipulations of the quantum system H-12]. This geomet- 
ric approach to quantum computing is argued to be more ro- 
bust against noise in the control parameters lfl3i[l4ll . Geomet- 
ric single- and two-qubit logic gates have been demonstrated 
in some systems such as nuclear spins Il5l4l7ll . trapped- 
ions 1 1_8] and superconducting qubits |T^, 20ll . For spin 1/2 
systems, theoretical calculations predict robustness against 
fluctuations in the static field and inhomogeneities in the mi- 
crowave field 12111 ; and this has also been explored experimen- 
tally with trapped ultracold neutrons M22I1 . 

Here, we demonstrate the implementation of a single-qubit 
geometric phase gate using adiabatic control of electron spins. 
We show that under current experimental conditions this leads 
to a much higher fidelity than achieved with simple dynamic 
phase gates. Interestingly, we find that the adiabatically ob- 



tained fidelity is comparable to that achieved by composite 
non-adiabatic pulses. These results are also verified by sim- 
ulations, which indicate that the fidelity of the adiabatic geo- 
metric phase gate remains high when the inhomogeneity in the 
microwave field strength becomes large, unlike for the non- 
adiabatic approach. 

After a cyclic evolution, a quantum system acquires a phase 
that depends on the geometric property of the evolution, the 
so-called Berry phase |8y. For an electron spin 1/2, this ge- 
ometric phase is determined by its trajectory on the Bloch 
sphere. Consider a spin initialised in the eigenstate |0) with 
respect to a static magnetic field along the z-axis, and with an 
additionally applied microwave field detuned from resonance 
(see Fig.Q]A). Then, slowly tuning the microwave frequency 
into resonance induces the eigenstate |0(i)} to adiabatically 
follow the effective magnetic field B{t) in the rotating frame, 
rotating it into the xy-plane. The phase of the microwave drive 
can be swept to rotate the eigenstate by some angle <fi. Fi- 
nally, the microwave field is detuned again, taking the eigen- 
state back to the z-axis. The geometric phase 7| ) acquired 
by the state |0) is given by the enclosed solid angle 9 of its 
trajectory on the Bloch sphere, 7| ) = 0/2 = <fi/2. By the 
same analysis, the geometric phase acquired for the |1) state 
is Til) = — 0/2> yielding a geometric phase for a general state 
of 7 = 7|0) - 7|i) = </>■ 

In addition to the geometric phase, the electron spin also 
acquires a dynamic phase during the evolution given by S = 
j-gp:B(t) ■ Sdt, where r is the total length of the control 
sequence and S is the expectation value of the spin operator, 
so that g/j,B(t) ■ S = E(t) is the time-dependent energy of the 
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eigenstate. To remove this dynamic phase from the final state, 
we introduce a tt phase shift of the microwave field exactly 
halfway through the control sequence. The dynamic phase 
accumulated during the second half of the control sequence is 
h = SI y^B(t) -Sdt = - j.gnB{t) ■ Sdt = -6 lf thus 
resulting in a vanishing final dynamic phase. 

There are two options for tuning the initially off-resonant 
microwave field into resonance: first, by adding an offset to 
the static field Bq field in ^-direction whose magnitude de- 
creases in time. Second, by tuning the frequency of the mi- 
crowave field. As the length of the adiabatic process is within 
a few microseconds, we employ the latter approach. 

For a time-dependent microwave frequency u>(t), the trans- 
formation of the Hamiltonian from the lab frame to a rotating 
frame with subsequent rotating wave approximation (RWA) 
can be performed in several ways. For instance, one can 
use the canonical rotating frame with constant frequency ui T f, 
where the spin resonance frequency ljq would be a natural 
choice for w r f . However, the resulting Hamiltonian then fea- 
tures important time-dependent oscillatory terms. Alterna- 
tively, we can choose a rotating frame which always tracks 
the frequency of the driving microwave field oj x { = uj(t). In 
this case, the transformed Hamiltonian after the RWA is given 
by (in the eigenbasis of S z ) 

'±(A + iA) Cle-tP 



H = h 



Qe lv -i(A-KA) 



(1) 



where fi, ip are the time-dependent amplitude and phase of 
the microwave field, A = u>(t) — ujq is the detuning, and 
A = uj(t) its time derivative. For our simulations this form 
of the transformation is used, ensuring that the measurement 
of the spin magnetization is carried out in the same reference 
frame that is used for the depiction of the trajectory on the 
Bloch sphere in Fig.[TjA. 

Fig. [TjB and C show an example for a microwave profile 
giving rise to an adiabatic geometric tt phase gate. The fast os- 
cillations at the beginning and the end of the control sequence 
in Fig.[TJC arise due to the finite detuning, getting slower as 
the microwave field approaches resonance (A — > 0), and dis- 
appearing entirely for td < t < td + t s . The variation of f2 x 
and Q, y during td < t < td + t s is due to the phase sweep of 
the microwave which causes the rotation of the spin magneti- 
zation in the xy-plane. 

Different choices for the microwave profile can be made, 
but vitally the time derivatives (A, Cl, tp) must be kept contin- 
uous (except at time t = td + t s /2) for achieving adiabatic 
evolution. However, continuity of the aforementioned param- 
eters is not sufficient for meeting the adiabaticity condition, 
which also requires that the control sequence r should take 
much longer than the spin precession, i.e. r 3> Aq , where 
Ao is the initial detuning. Experiment and simulation both 
confirm that with all the other parameters fixed, the fidelity of 
the adiabatic phase gate then increases for larger values of r. 

We now simulate such an adiabatic tt phase gate and 
compare its performance with phase gates based on dy- 
namic pulses. The adiabatic implementation is based on 
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FIG. 1: A: Evolution of the spin eigenstate |0) (red trace Ci) repre- 
sented by the spin vector S and effective magnetic field arising from 
the applied microwave field (blue trace C2). B: The detuning A, am- 
plitude il, and phase tp of the microwave field during an adiabatic n 
phase gate. The microwave frequency is tuned from off-resonance to 
resonance between time and td and then kept resonant for a period 
of t 3 while the phase <f> of the microwave is swept. After time td + t s 
the microwave is detuned again. The effective field changes sign af- 
ter the 7r phase shift at td + t 3 /2. C: The x (blue) and y (green) 
components of the microwave field applied to the electron spin for a 
7T phase gate. 



the pulse sequence shown in Fig. Q] C, whereas the se- 
quence for applying a dynamic phase <j> is {W) x — (4>) y — 
(w)-x- The corresponding BB1 composite pulse sequence 
is built by replacing each single pulse with a compos- 
ite pulse (0) lpo (ir) tpi+lpo (2ir) ip2+ipo ('K) !pi+!po , where 6 and 
<Po G {x, y, —x} are the rotation angle and phase of 
the corresponding single pulse, respectively, and ipi — 
arccos(— 9/Att), tp2 = 3<pi. It is also possible to construct a 
non-adiabatic geometric phase gate by applying two tt pulses 
of different phase successively (tt) IPi — (tt) V2 , and the BB1 
version of this sequence can be built accordingly. The geo- 
metric phase acquired by this operation is determined by the 
phase difference between the two tt pulses: 7 = 2(<p2 — <pi). 

The simulated phase errors of the three different phase gates 
are plotted as a function of the deviation from a given mi- 
crowave amplitude (Fig.|2]A). The maximum relative devia- 
tion Afio from its center value fio is set to ±10%, a common 
value for typical ESR spectrometers. The results show that 
the phase error in the BB1 dynamic phase gate is much larger 
than for both geometric phase gates. The relative phase er- 
ror has a third order dependence on the error in microwave 
amplitude for the BB1 dynamic phase gate, i.e. e r dyn ~ 
O [( AQq / r2o) 3 ] , whereas it has a sixth order dependence for 
the BB1 geometric phase gate, i.e. e r .g CO ~ O [(Ailo/^o) 6 ] ■ 
The relative phase error in the adiabatic geometric phase gate 
is comparable with the BB 1 geometric phase gate in this range 
of microwave inhomogeneity. However, we can see that if the 
microwave field strength inhomogeneity becomes larger than 
10%, which is the case for most coplanar cavities, the per- 
formance of the adiabatic phase gate exceeds that of the BB 1 
geometric phase gate. Small fluctuations of the adiabatic gate 
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FIG. 2: A: Simulated relative phase errors e r = 'y/n — 1 for tt phase 
gates as functions of the error in microwave amplitude Afio/^o, 
where 7 is the acquired phase in the simulation. B: Infidelities I 
(see text) of the operations as functions of ASlo /fio- All simulations 
employ a single spin 1 /2 and do not include spin packet effects. The 
parameters are fio = 14 MHz, initial detuning Ao = 6 MHz, td = 
2 [xs, t a = 4 /is. Note that since there is no microwave noise, the 
BB1 errors are independent of the BB1 pulse duration. 



are due to imperfect adiabaticity of the operation, decreasing 
both for a stronger microwave field or a slower passage. 

The quality of a phase gate can also be characterised by the 
infidelity / of the operation, defined as 



1=1- 



iTrlC/CV 1 ] 



(2) 



where Uq and U are the operators for an idealized and simu- 
lated tt phase gate, respectively. Fig.|2]B shows the infideli- 
ties of the three phase gates as functions of AOo/fio- F° r the 
two phase gates using BB1 composite pulses, the infidelity 
increases as O [(Afi /^o) 6 ] ■ The infidelity of the BB1 dy- 
namic phase gate is much smaller than its phase error, because 
in the measurement of infidelity only the diagonal elements of 
U are considered while the phase error e r is more directly re- 
lated to the off-diagonal elements of U. The oscillations in 
the infidelity of the adiabatic phase gate are due to the finite 
non-adiabaticity of the operation, but only have a minor effect 
on the geometric phase. Once again, the infidelity of the BB1 
geometric gate is comparable with the adiabatic variant within 
±10% microwave inhomogeneity, yet the adiabatic phase gate 
is more robust when the inhomogeneity becomes larger. 

For the experiments, we used a sample with narrow ESR 
linewidth (P donors in high-purity 28 Si crystal at 8K) in or- 
der to ensure that all the spins are within the bandwidth of the 
adiabatic control sequence. The X-band microwave signal is 
generated at a constant frequency, which is then modulated 
by the I/Q signals from an arbitrary waveform generator to 
create the required microwave field, such as the one shown 
in Fig. [TJC. The complete sequence for measuring the geo- 
metric phase gate consists of an initial (dynamic) tt/2 pulse 
that creates the spin coherence in the xy-plane, the adiabatic 
control sequence, and a (dynamic) tt pulse that refocuses the 
random fluctuations of the environment (Fig. [3] A inset). The 
spin echo is detected and its phase is determined by quadra- 
ture detection, from which the phase acquired by the electron 
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FIG. 3: A: Measured geometric phase 7 (black) and error e 7 (red) as 
functions of the angle <f> that the spin rotates in the xy plane during 
the adiabatic sequence. Inset: Pulse sequence for measuring the ge- 
ometric phase gate. B: Intensities of the measured spin echo signal 
normalized to corresponding Hahn echo intensities. C: Phase error 
and D: Echo intensity for adiabatic phase gates of different t s . In 
these experiments td = 2 fis, maximum of the microwave amplitude 
fio = 0.64 MHz, and initial detuning A = 6 MHz. 



spin during the adiabatic phase gate can be deduced. 

Fig. [3] A shows the phase of the electron spin 7 measured 
after an adiabatic phase gate that is designed to apply a geo- 
metric phase to the spin, and the corresponding error defined 
as e 7 =7 — 0. The experimental data follows the theoreti- 
cal relation 7 = very well over a broad range of 0, from 
to 207T, which verifies that we have successfully implemented 
the adiabatic geometric phase gate to the electron spins. The 
intensity of the spin echo is also plotted against (Fig.[3]B) as 
a measurement of the performance of the phase gate. The echo 
intensities are normalized to a Hahn echo for the same time 
delay. The fact that the echo intensity at = is less than 1 
indicates that the adiabatic process is not perfect, and implies 
only partial adiabatic following of the whole spin ensemble. 
This is partly due to the off -resonance error of the spins, how- 
ever, since the ESR linewidth of the spins is narrower than the 
bandwidth of the adiabatic control sequence, the failure of adi- 
abatic following is more generally due to the non-adiabaticity 
of the phase gate operation. This also explains why the echo 
intensity decreases for greater 0: for a fixed duration of the 
adiabatic sequence; a greater implies a faster phase varia- 
tion in the xy plane during td < t < td + t s , hence a less 
adiabatic operation. In addition, because of the off -resonance 
error and inhomogeneities in the microwave field, different 
spin packets do not follow exactly the same path during the 
adiabatic operation, and by the end of the evolution they will 
exhibit a spread in the final phase. While the geometric phase 
7 is a measurement of the mean phase of all the spins, the spin 
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FIG. 4: A: Measured phase error and B: echo intensity for the adi- 
abatic geometric phase gate, and dynamic phase gates using single 
microwave and BB1 composite pulses. All intensities are normal- 
ized to unity to enable a more convenient comparison. The parame- 
ters used for the adiabatic sequence of the experiment are the same as 
in Fig. [3] C and D: Simulation using the parameters On = 1.4 MHz, 
An = 6 MHz, t d = 2 fis, t a = 4 /is. 



echo intensity reflects the variance of phases of between dif- 
ferent spins. The reduction of the echo intensity can thus be 
attributed to the loss of phase coherence of the spins. 

The effect of varying t s , which determines the sweep rate 
of <j> in the xy plane, is shown in Fig. [3] C, D. In accordance 
with the adiabatic condition, the phases measured with shorter 
t s contain larger errors. Furthermore, the echo intensity is re- 
duced for shorter t s implying higher non-adiabaticity. The 
dips in the echo intensity traces for t s = 3 fis and 4 fis are 
attributed to the non-adiabaticity of the phase gate rather than 
noise from the spectrometer, as we have observed similar fea- 
tures in simulations where the only imperfection introduced is 
the microwave field inhomogeneity. 

We proceed by comparing the adiabatic phase gate to non- 
adiabatic phase gates based on single microwave and BB1 
composite pulses. The gates are studied in the range of [0, An] 
since both the BB 1 pulse operation and the non-adiabatic ge- 
ometric gate have a natural phase limit of 4ir. The pulse se- 
quences for the non-adiabatic gates are as described above for 
the simulations. Limited by the output level of the signal gen- 
erator, the range of the linear amplification of the spectrom- 
eter, and our solid-state amplifier, the maximum amplitude 
of the microwave field we can apply is about 0.25 G, corre- 
sponding to the length of a it pulse = 700 ns. In this 
case, the BB1 pulse sequence for a dynamic it phase gate is 
14 x = 9.8 fis, therefore we choose t s = 6 fis and td = 2 fis 
for the adiabatic sequence so that its total length t = 10 /is is 
comparable with the BB1 pulse sequence. 



The measured phase error and echo intensity shown in 
Fig- SI A and B, respectively, demonstrate that the adiabatic 
phase gate outperforms its dynamic counterpart using single 
microwave pulses, while its performance is comparable to a 
dynamic gate with BB1 composite pulses. The comparison 
between experimental data and simulations (Fig. |4] C and D) 
also verifies that the adiabatic geometric phase gate is more 
robust under B\ inhomogeneity. We have performed a similar 
comparison of the adiabatic and the non-adiabatic geometric 
phase gate, which gives the same qualitative results, i.e. BB1 
pulses need to be employed for the non-adiabatic implemen- 
tation to obtain comparable performance. For the simulations, 
we assumed 10% inhomogeneity in B\ and no off -resonance 
effect. In Fig. |4]D the simulation reproduces the main fea- 
tures in the echo intensity of the dynamic phase gate using 
single pulses (red trace), which implies that the reduction in 
the echo intensity is essentially due to the B\ inhomogeneity. 

A previous theoretical study has suggested that a slow adia- 
batic process is more exposed to environmental decoherence, 
mitigating its advantage over non-adiabatic operations 123 1. 
However, in our experiment the gate times for the BB 1 and 
adiabatic phase gate are similar, and the fidelity of the adi- 
abatic geometric phase gate is still limited by other imper- 
fections of the equipment such as the small amplitude of the 
microwave field and phase imprecision. 

In conclusion, we have demonstrated single-qubit geomet- 
ric phase gates using adiabatic control of electron spins. Ex- 
periments and simulations showed that the adiabatic geomet- 
ric phase gate is remarkably robust against inhomogeneities 
in the microwave field. The performance of the adiabatic ge- 
ometric phase gate should improve further with a larger mi- 
crowave field amplitude, for example achieved using a higher 
Q-value resonator 12411 . For the current experimental setup its 
performance is comparable to the geometric phase gate using 
composite non-adiabatic pulses such as BB1 pulses. How- 
ever, the adiabatic phase gate could be advantageous given a 
more inhomogeneous microwave field (such as may arise in 
coplanar resonators 12511 ), or at higher microwave amplitudes. 
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